Spin-Correlation Coefficients and Phase-Shift Analysis for p+^He Elastic Scattering 
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Angular distributions for the target spin-dependent observables Aoy, A^x, and A^j, have been 
measured using polarized proton beams at several energies between 2 and 6 MeV and a spin-exchange 
optical pumping polarized '^He target. These measurements have been included in a global phase- 
shift analysis following that of George and Knutson, who reported two best-fit phase-shift solutions 
to the previous global p-|-^He elastic scattering database below 12 MeV. These new measurements, 
along with measurements of cross-section and beam-analyzing power made over a similar energy 
range by Fisher et al, allowed a single, unique solution to be obtained. The new measurements 
and phase-shifts are compared with theoretical calculations using realistic nucleon-nucleon potential 
models. 

PACS numbers: 21.45.-v, 21.30.-x, 24.70.+S, 25.40.Cm 



I. INTRODUCTION 

Beginning from the picture that atomic nuclei are 
composed of interacting nucleons, ab initio calculations 
of light nuclear systems are based on realistic nucleon- 
nucleon potential models, which have been adjusted to 
reproduce two-nucleon (NN) scattering and bound-state 
data accurately [l|. This effort has included calculations 
of low-energy scattering observables for the three-nucleon 
(3N) and four-nucleon (4N) systems. The latter calcula- 
tions are especially significant because the 4N system is 
the lightest to exhibit thresholds and resonances [Ij, so 
that its correct description is an important milestone for 
this approach. 

The comparison of ab initio calculations with nucleon- 
deuteron scattering measurements reveals general agree- 
ment for the cross-section and tensor analyzing powers, 
but significant underprediction of the beam and target 
analyzing powers [3]. A similar "Ay Puzzle" has been 
reported for p+'^He elastic scattering [i^ . All NN models 
and theoretical methods yield this disagreement, which 
is not resolved by including the 3N force necessary to 
reproduce the 3N and 4N binding energies 

The study of this discrepancy may benefit from the 
more thorough comparison between theory and experi- 
ment made possible by a set of experimental phase-shifts 
and mixing parameters. For example. Ay is known to 
be particularly sensitive to the splitting between triplet 
P-wave phase-shifts |4J] . A wealth of experimental data 
exists for p-|-^He elastic scattering below 12 MeV pro- 
ton energy, with the niost recent phase-shift analysis by 
George and Knutson jQ] performed on a database of over 
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1000 data points. That analysis, however, was unable 
to constrain a unique set of parameters, and instead ob- 
tained two solutions which fit the data equally well. The 
difference between the two solutions was largest for spin- 
correlation coefficients below 4 MeV, where no such data 
exist. With the aim of resolving the phase-shift ambigu- 
ity, we have used a polarized '^He target Q to measure 
angular distributions of the spin-correlation coefficients 
Axx and Ayy at proton energies between 2 and 6 MeV 
and included those new data, along with those of Fisher 
et al. in a new global phase-shift analysis following 
that of George and Knutson. 



II. EXPERIMENT 
A. Polarized Target 

Polarized and unpolarized beams from the Triangle 
Universities Nuclear Laboratory (TUNL) tandem accel- 
erator were directed by an analyzing magnet to a 62 cm 
diameter scattering chamber. The polarized '^He target 
was specifically designed for low-energy charged-particle 
scattering experiments and has been described previously 
in detail [7|. In contrast to previous targets used for the 
same purpose [sJ-fT^] which operated by metastability- 
exchange optical pumping (MEOP), the present tar- 
get used spin-exchange optical pumping (SEOP). This 
method has previously been used in polarized ^He tar- 
gets for electron [l^ and gamma-ray [l^ scattering ex- 
periments. 

The primary advantage of SEOP over MEOP is greater 
target thickness, since it polarizes ^He at a pressure of 
about 8 bar instead of several mbar. The need to mini- 
mize energy loss for the incident and scattered particles 
for the low-energy application, however, required the use 
of thin windows to contain the gas. Since such windows 
cannot withstand the full 8 bar '^He pressure, we opti- 
cally pumped '^He in a system separate from the scatter- 
ing target cell. The target cell was then batch-filled with 
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polarized gas to a pressure of approximately 1 bar. 

The target cell was a 5.1 cm Pyrex sphere with open- 
ings along the equator for the incident and scattered par- 
ticles. These were covered with 7.6 ^m Kapton foil af- 
fixed with Torr Seal epoxy The cell was housed in 
a compact sine-theta coil to provide a uniform 0.7 mT 
magnetic holding field. An NMR coil pressed against the 
rear of the cell was used to measure the target polariza- 
tion, as discussed below. 

The ^He polarization was produced by SEOP using Rb 
as the intermediate alkali metal. A 60W fiber-coupled 
diode laser system tuned to the 795nm Rb Dl absorp- 
tion line provided, with appropriate optics, the circularly 
polarized light for the optical pumping. Two modifica- 
tions to our original polarizer [7] were attempted. In 
agreement with results reported by others |17l] . the use 
of "mixed metal" optical pumping cells containing both 
Rb and K was found to decrease the "spin- up" time [l^ ■ 
In this work, the typical time required to reach satura- 
tion polarization was about 12 h with mixed-metal cells, 
compared to about 24 h using only Rb. 

The other modification was the use of frequency- 
narrowed laser light for optical pumping. Following 
the work of [l9j . an external Littrow cavity was con- 
structed and used to reduce the output width of 50 W 
Quintessence (20| diode bar array from 2 nm to 0.3 nm 
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The narrowed output power was about 30 W. Un- 
fortunately, the ^He polarization produced with this laser 
system was not consistently higher than that produced 
with the 60-80 W broadband system, perhaps because 
more light was absorbed from the latter. The majority 
of the present data was therefore taken with the broad- 
band laser. 



B. Scattering Measurements 

The experimental arrangement is illustrated in Fig. [1] 
Beam current on the target cell was limited to 50 uA 
to minimize damage to the Kapton foils. Failure of the 
epoxy, especially that sealing the beam-exit foil, caused 
cells to leak after a few days. Each such failure required 
a cell change and recalibration of the NMR signal (see 
below). 

Measurements were made at five proton energies be- 
low 6 MeV which overlap the e nerg ies of both Fisher et 
al. U] and Alley and Knutson [Tj]. Analyzing magnet 
settings which determined the beam energy were selected 
according to the calibration of Ref. ^3] . Beam energies 
were adjusted to offset e nerg y loss in the foils and gas, 
as modeled with TRIM 
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data taken with different thickness entrance foils were 
slightly different, and the error- weighted average value 
was adopted. The uncertainty in TRIM stopping pow- 
ers for materials used was estimated by comparison with 
experimental stopping powers (23j . and ranged from 3- 
10%. An uncertainty of 10% was assigned to cases where 
no data were present in the relevant energy range. 



The beam and target polarizations were reversed fre- 
quently during data-taking. The beam polarization was 
reversed at either 1 or 10 Hz in the sequence "udduduud" , 
where "u" means "spin-up" and "d" means "spin-down". 
The target polarization was reversed less frequently, since 
a few seconds were required to reverse the target's mag- 
netic field. Polarized target data were collected for inter- 
vals of 2.5 m in each spin orientation, with NMR mea- 
surements of the magnitude of the polarization made im- 
mediately before and after the orientation was reversed. 
The target polarization decayed with a 2-3 h time con- 
stant, so this process was stopped when the gas was 
judged to be too depolarized, generally after about 1 h. 
At that time the gas was exhausted from the target cell, 
which was then flushed with research-grade N2 and re- 
filled with a new batch of polarized gas. The recovered 
depolarized ^He gas was circulated through a LN2-cooled 
trap to remove impurities before repolarization. 



Scattered particles emerging from the target were de- 
tected by four pairs of Si detectors which could be rotated 
to the desired angle. Available angles were restricted by 
the windows in the sine-theta coil's mu-metal shield to 
20° increments between 30 and 150°. The shield could 
be moved axially so that "intermediate" angles offset by 
10° were also available. The detectors were each placed 
in an Al holder behind two brass collimators spaced 5.08 
cm apart in an Al "snout" , which restricted the range 
of scattering angles visible to the detector to 1.5°. De- 
tectors were as close as possible to the target without 
the 30° snout's touching the sine-theta coil, so that the 
distance from the center of the target cell to the front 
collimators was about 10.2 cm. Beam current on target 
was measured by a electrostatically suppressed Faraday 
cup located about 0.5 m behind the target cell. Charge 
went to ground through a current integrator to measure 
the relative number of beam particles in each spin state. 



An example detector spectrum is shown in Fig. [5] The 
peaks corresponding to protons elastically scattered from 
■^He and small amounts of N2 were well separated. The 
latter was required by the optical pumping process 0. 
Occasionally *He was added to the cell to measure the 
beam polarization, and in some cases that peak over- 
lapped the ■^He peak. In such cases the peaks were fit 
with skewed Gaussians. A pulser was added to each 
spectrum to measure electronic dead time, which was 
typically less than 1%. 



The observables were extracted from peak yields in 
left (L) and right (R) detectors using an extension of 
the geometrical mean method (23 | for analyzing powers 
to include polarized beam and target. With the spins 
aligned vertically along the ±y-axis, the following cross- 
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FIG. 1: Diagram of the experimental arrangement. The polarimeter chamber at left was used with the detectors either 
horizontal, as shown, or vertical. The polarimeter cell was removed from the beam during data taking to allow beam to reach 
the target chamber. The coordinate system used to define the scattering observables is also shown. 
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FIG. 2: ='He(p,p)^He spectrum taken at 5.54 MeV and 90° 



for example L refers to the number of particles scat- 
tered into the left detector while the beam was spin "up" 
and the target was spin "down" . The detector yields 
were normalized to the current integration and target 
pressure, which was known to about 2%, for each spin 
state. The polarization of the beam and target are given 
by p and ^ respectively. The observables are therefore 
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Similarly, when the beam and target spins are aligned 
horizontally along the x-axis. 
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If either the beam or target is unpolarized, only one 
analyzing power will be non-zero, and its expression re- 
duces to the usual cross-ratio for analyzing powers. 



ratios were formed: 



Xi 



Xo = 
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The same is true for the scattering of spin- 1/2 protons 
from spin-0 a-particles used for beam polarimetry. 



C. Beam Polarimetry 



Ln+LU\ fRn +R 1 + PyPi ^yy 



^LU+LitJ \Rn+RitJ ~ l-PyP^Ayy 



where the arrows indicate the beam and target spin state. 



Proton beam leaving an atomic beam polarized ion 
25| passed through a calibrated Wien filter at 



source 



the ion source to orient the spin quantization axis of the 
beam in the desired direction at the scattering chambers. 
The magnitude of the beam polarization was measured 
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periodically using p+ He elastic scattering in either the 
target cell or in a separate cell in a polarimeter cham- 
ber installed upstream of the target chamber. Detectors 
in the latter could be mounted at 110° in the horizontal 
or vertical scattering planes, so that either polarization 
component could be measured. The cell in the polarime- 
ter chamber could be moved to insert it periodically into 
the beam for polarization measurements. 

The p-|-*He asymmetries were divided by the analyz- 
ing power Ay to obtain the beam polarization. Pub- 
lished phase shifts [l^ were used in a spini-on-spin-0 
phase-shift code to calculate the analyzing powers for 
the energies at the center of the cell as determined from 
TRIM. The uncertainty in the resulting analyzing power 
was typically 2%. 

For more than half of the spin-correlation data, how- 
ever, the beam polarization was unstable, so that periodic 
monitoring did not necessarily determine the average po- 
larization. Therefore, the beam polarization for all Aj,j, 
measurements was determined by normalizing our rela- 
tive Aj^o measurements to published values flSl, l22l| . Each 
point in a relative Aj,o angle set was divided by a value 
linearly interpolated from those previous measurements 
at the same energy, and the polarization was taken to 
be the average of these ratios. An uncertainty of 0.02 
was assigned to the polarization and added in quadrature 
with statistical uncertainties. No published data were 
available at 2.7 MeV, so smooth curves were fit through 
existing distributions at each angle vs. energy and eval- 
uated at 2.7 MeV. The normalization then proceeded as 
above. 

The procedure was extended to about one-third of 
the measurements by "tipping" the spin 20° out of 
the plane with the Wien filter and applying the above 
analysis to the y-component. The x-component of the 
beam polarization was obtained by multiplying the y- 
component by the ratio of the two components. An 8% 
normalization uncertainty was applied to those angular 
distributions to account for the estimated 1.5° uncer- 
tainty in the relative azimuthal orientation of the scatter- 
ing plane and the polarized beams quantization axis. The 
remaining Kxx measurements with stable beam polar- 
ization relied on polarimeter measurements as described 
above. 



D. Target Polarimetry 

As discussed in detail by Katabuchi et al. p7|, the 
target polarization was monitored using pulsed NMR. 
Briefly, an RF pulse at the Larmour frequency set by the 
sine-theta coil magnetic field was sent through a small 
coil pressed against the rear of the target cell. The re- 
sulting collective precession of "^He spins about that mag- 
netic field induced a signal proportional both to the '^He 
polarization and pressure. This voltage was then divided 
by the cell pressure to give a relative measure of polar- 
ization. 



These relative NMR data for each target cell were cal- 
ibrated against separate "'He -I- ■^He Ay measurements at 
an energy and angle where Ay= -1. The resulting scatter- 
ing asymmetries for an incident "^He beam, given by Eq. 
[5l were therefore direct measurements of the "^He target 
polarization. This calibration method was motivated by 
the prediction of Plattner and Bacher [27j of an Ay = -1 
extremum near 15.33 MeV '^He lab energy and 47° '^He 
lab scattering angle. Their prediction for its location was 
only approximate. Thus, relative measurements of Ay in 
'*He-|-'^He elastic scattering as a function of angle and en- 
ergy near the predicted extremum were made to define 
the local minimum. We determined 6„im to be 46.64 ± 
0.22°, and the measurements of Ay vs energy agree with 
the prediction of Ref. [28^ of a very broad minimum. 

The simultaneous NMR-scattering calibrations for var- 
ious target cells were made at 45°, and Eq between 15.44 
and 15.82 MeV. Since this was not exactly the minimum 
point for Ay, a target polarization value was assigned 
from the relative measurements with an uncertainty of 
3%, normalized so that the minimum Ay was equal to -1. 



E. Steering Effect 

The target magnetic field, though small, steered the 
incoming and scattered protons slightly through the 
Lorentz force, as sketched in Fig. [3l When the B-field 
was reversed to reverse the target spin, the particles were 
steered oppositely, and the relationship between detector 
yields reversed. 

First, the incoming beam deflection moved the scat- 
tering center closer to one detector and further from the 
other, thus changing their relative solid angles. Second, 
the scattered protons were deflected to emerge from the 
target cell at slightly different angles than those at which 
they were actually scattered. Thus, for a given orienta- 
tion of the magnetic field, the actual scattering angle 
of particles reaching one detector of a left/right pair was 
more forward than the detector angle, while that of those 
reaching the other detector was more backward. This 
difference in scattering angle produced an instrumental 
asymmetry through the angular dependence of the dif- 
ferential cross section. These two steering effects, unlike 
other systematic effects [s^l , produced systematic instru- 
mental asymmetries which were not cancelled by revers- 
ing the target spin. 

These effects were largest at our lowest bombarding 
energies. Fig. 2] shows observables extracted from data 
taken while the target was unpolarized, when the asym- 
metries corresponding to the target analyzing power and 
spin-correlation coefficients should be zero. While this 
was true for Aqj, measured with the sine-theta coil's 
magnetic field turned off, non-zero asymmetries were ob- 
tained with the field on. The asymmetries for A^y and 
Axx^ on the other hand, were consistent with zero even 
with the field on. 

The result of a simple calculation of these two effects 
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of the observable at forward angles at the lowest proton 
energies. 



Beam 



Results for Observables 



Right Detector V 

FIG. 3: Steering of the incoming proton beam and scattered 
particles by the sine-theta coil magnetic field. The target cell 
and detectors are shown from above. The beam is incident 
from the left on the gas cell along the dotted line, but deflected 
as shown by the magnetic field, which is oriented out of the 
page. Similarly, the scattered particles travel along curved 
paths to the detectors. The figure is not to scale, and the size 
of the effect is exaggerated for clarity. 
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FIG. 4: (Color online) Instrumental target asymmetries mea- 
sured at 2.25 MeV with the sine-theta coil magnetic field both 
on and off. A simple calculation of the effect of magnetic 
steering is also shown. 



at 2.25 MeV is also shown in Fig. SI Though the largest 
deflection angle calculated was 0.1°, the resulting asym- 
metries are large enough to interfere with polarized target 
measurements. The calculation neglects the finite size of 
the beam and target, and simply determines the energy 
loss, modeled in TRIM, and magnetic steering of incom- 
ing and scattered protons in small steps as they proceed 
through the magnetic field. The calculation reproduces 
the general size and forward-angle trend of the effect, but 
with insufficient detail to be used to correct the data. 

Instead, actual measurements of these instrumental 
asymmetries, obtained both from direct measurements 
with unpolarized ^He and by extrapolating polarized tar- 
get asymmetries to zero polarization, were subtracted 
from polarized target Aq^ data. The corrections obtained 
in this way were often large, being several times the size 



The present measurements of Aqj,, A^j,, and K^x are 
shown in Figs. [5]l7l Each plot also includes a curve 
calculated from the new best-fit effective range param- 
eters of the global phase-shift analysis (PSA) discussed 
below, as well as previous measurements where available. 
The overall agreement with the previous measurements 
is good, with the most forward-angle A^j, points of Al- 
ley and Knutson [l^ at 5.54 MeV being the only excep- 
tion, and the present results have smaller error-bars. The 
new measurements are well-fit by the phase-shift analy- 
sis, except for the two most forward-angle Aqj, points at 
3.15 and 4.02 MeV. The forward-angle points required 
the largest correction for magnetic steering, so the dis- 
agreement with the PSA may indicate that the correction 
applied to those points was not sufficiently accurate. 



III. PHASE-SHIFT ANALYSIS 

A phase-shift analysis of the global p+'^He elastic scat- 
tering database below Ep= 12 MeV was performed fol- 
lowing the earlier work of George and Knutson [6] , with 
the addition of about 300 new data points, including the 
^ and Aj,o measurements of Fisher et al. and the 
present Aoy, A^j^, and Plxx measurements. These addi- 
tional data all fell between 1.0 and 5.54 MeV. The search 
routine was the same as that used in the previous analysis 
and was provided by George [sl] . 

The program calculated scattering observables as func- 
tions of scattering matrix elements, which were in turn 
parameterized using phase-shifts and mixing parameters 
according to the Blatt-Biedenharn convention [s^. The 
phase-shifts and mixing parameters used were "'^Sq, '^Si, 
iPo, 3P2, 3Pi, 3Po, ID2, £ (1-), e (1+), and e (2"), as well 
as consolidated and '^Yj triplet phase-shifts. The en- 
ergy dependence of the phase-shifts and mixing parame- 
ters was described by the first three terms in a modified 
effective-range expansion. These 36 effective range pa- 
rameters were adjusted to minimize with respect to 
the experimental database using the MINUIT package 
[33 |. As described in Ref. j6|], the database was broken 
into groups of measurements thought to have common 
normalizations. These 21 normalization factors were an- 
alytically adjusted at each step of the parameter search 
to further minimize x^. 

Initial parameter searches resulted in multiple solu- 
tions which were discontinuous in one or more of the 
^D2, '^Dj, and '^F^- phase-shifts. The discontinuity, dis- 
cussed also by Alley [IHl, occurred when the phase-shift 
crossed zero and was deemed unphysical. The number of 
such solutions was reduced by fixing the small '^Fj phase- 
shift at the values obtained using the database of Alley 
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FIG. 5: (Color online) Present measurements of Ao^ at five 
proton beam energies, together with the global PSA fit. Mea- 
surements of Szaloky et al. [3l[| and Alley and Knutson [T^ 
as well as theoretical calculations of Deltuva [s^l using the 
CD-Bonn potential ^SQ] are also shown. 



FIG. 6: (Color online) Present measurements of Ayy at five 
proton beam energies, together with other results as men- 
tioned for Fig. 8. 



and Knutson [T^. The best- fit when this parameter 
was not searched increased by only 0.2%. All but one of 
the remaining solutions, which spanned a range of about 
10% in x^, were rejected by demanding that the phase- 
shifts be continuous in the energy range covered by the 
database. The remaining solution, which has the lowest 

and is adopted as the present global result, yields a 
small positive scattering length for ■^Dj , indicating that a 
discontinuity must occur in that phase at an energy be- 
low the lower end of the database (100 keV). The best-fit 
effective range parameters and associated statistical un- 
certainties are given in Table HI 

The global solution had a x'^-per-datum of about 2 
for the data added in this analysis. This could be im- 
proved to between 1.3 and 1.5 if points whose individual 

contributions exceeded 10 were rejected. About half 



TABLE I: Best-fit effective-range parameters with statistical 
uncertainties. 



Phase 


aOxlQ-^ 


alxlO"' 


a2 


^So 


-9.0±0.4 


7.9±0.6 


1.0±0.2 


^Si 


-11.06±0.17 


7.5±0.3 


-0.09±0.08 


^Pi 


5.44±0.15 


-1.7±0.3 


4.05±0.13 


^P2 


2.128±0.012 


1.591±0.018 


0.356±0.009 


^Pl 


1.63±0.02 


2.06±0.04 


0.592±0.018 


^Po 


8.8±0.3 


-1.5±0.3 


2.40±0.09 


'B, 


-14±4 


15±4 


-8.3±1.3 


^D, 


-0.06±0.19 


-3.4±0.9 


7.2±0.4 


^F, 


4.19 


38.3 


-6.29 


e(l+) 


-5±5 


34±6 


-10.7±1.7 


e (1") 


-420±11 


230±10 


-41±2 


e (2") 


17±7 


8±6 


-5.9±1.4 
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TABLE II: Global phase-shift analysis results at 2.25 and 3.15 
MeV 





2.25 MeV 


3.15 MeV 


Phase 


Present 


CD-Bonn 


Present 


CD-Bonn 




-39.1±1.7 


-39.6 


-48.7±0.9 


-49.3 




-34.5±0.7 


-34.8 


-42.90±0.09 


-42.9 


'Pi 


8±2 


10.6 


13.4±0.4 


14.9 




5±6 


7.9 


9.7±0.8 


12.3 


'Pi 


17±4 


16.9 


27.0±1.9 


26.1 




16.5±0.7 


16.0 


27.7±1.2 


25.8 




-10±20 


-8.9 


-12.2±1.7 


-8.3 
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FIG. 7: (Color online) Present measurements of A^x at three 
proton beam energies, together with other results as men- 
tioned for Fig. 8. 



TABLE III: Global phase-shift analysis results at 4.00 and 
5.55 MeV 





4.00 MeV 


5.55 MeV 


Phase 


Present 


CD-Bonn 


Present 


CD-Bonn 


^So 


-56.3±0.6 


-56.8 


-67.8±0.9 


-67.1 


''Si 


-49.3±0.5 


-49.7 


-58.6±0.3 


-59.2 


'Pi 


17.3±1.6 


18.2 


21.2±1.7 


22.5 


'■'Po 


14.1±0.9 


16.6 


21.3±0.7 


23.9 


■'Pi 


34.9±0.3 


33.9 


45.2±0.5 


43.0 


'P2 


37.6±0.6 


34.9 


51.5±0.5 


47.0 




-13±2 


-9.0 


-14±2 


-9.6 



of these 13 out of about 300 new points seemed simply 
to be random outliers, while the others seemed to be as- 
sociated with systematic problems. These included some 
forward angle Aq^ points which had been corrected for 
magnetic steering. Another apparent systematic problem 
was found for Aj,o data at 1.60 MeV for which the four- 
most-forward angle points disagreed with the phase-shift 
analysis. The effect of their removal on the phase-shifts 
was generally negligible, and in all cases within the range 
of systematic error indicated by the single-energy analy- 
ses described below. 

In order to gauge the effects of systematic errors, 
single-energy analyses were performed at energies where 
new spin-correlation and new or existing cross-section 
measurements were available, i.e. at nominal proton en- 
ergies of 2.25, 3.13, 4.00, and 5.54 MeV. AU measure- 
ments within 100 keV of the nominal energies were in- 
cluded. The same method was used for these single- 
energy fits as for the energy-dependent work, except 
that the phase-shifts were searched directly, instead of 
through the effective range parameters. 

The present phase-shift results, both single-energy and 
global, are shown in Fig. [8l along with those of Refs. [36| 
and 6] . The addition of the new data removes the S-wave 
ambiguity in the latter results without qualitatively mod- 
ifying the behavior of the other parameters, such as the 
resonant P-wave behavior associated with excited states 
of^Li. The new, low energy data also seem to introduce 
some tension with previous, higher energy data, as in- 



dicated by differences between the present global results 
and those of Ref. ^ for ^Sq and ^Pg. 

Results for selected phase-shifts, in degrees, are tab- 
ulated in Tables HHIIIII at the nominal energies of the 
single-energy analyses. The quoted uncertainties include 
statistical uncertainties and systematic sources including 
differences between the single-energy and global analyses 
and variation of the parameters when outliers are ex- 
cluded. Theoretical calculations, described below, using 
the CD-Bonn realistic NN potential are also shown for 
comparison. 



IV. COMPARISON WITH THEORETICAL 
CALCULATIONS 

We first compare our new experimental results with 
recent ab initio momentum-space calculations from Del- 
tuva and Fonseca [37[ which rigorously include the 



Coloumb interaction and use a variety of 2N potentials. 
For simplicity, only their results obtained using the CD- 
Bonn potential [38] are shown in Figures [5][71 but the 
results of the other realistic 2N potentials considered, 
AV18 [39] and N3L0 are similar. The Aq^ calcu- 
lations consistently underpredict the new results by 10- 
20% at the maximum. This is similar to, though smaller 
than, the previously established 40% underprediction of 
Ayo by several realistic potentials and theoretical meth- 
ods [I, [13. 

The theoretical results for the spin-correlation coeffi- 
cients at Ep = 2.77 MeV and above agree with the present 
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FIG. 8: Phase-shift results, in degrees, as functions of proton lab energy in MeV. Both the global and single-energy results of 
this analysis are shown, as well as the previous results of Alley and Knutson 36] and George and Knutson [^. 



results at backward angles but are too small by about 
0.02 at forward angles. The disagreement for backward 
angles between the theoretical results and the present 
phase-shift analysis at 2.28 MeV may result from the 



lack of back-angle data points, especially for K^x- 

For the phase-shifts, the theoretical S-waves generally 
agree well with the present results, while the theoreti- 
cal ^Pi and ■^Pq phase-shifts are larger. The theoretical 
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•^Pi and phase-shifts, as weh as the e(l~) mixing 
parameter, are consistently smaher than our present re- 
sults. The splitting between the triplet P- waves is also 
underpredicted, as shown in Fig. [3J where S — '^P2 - ( 
•^Pi -I- •^Po)/2. This is interesting in light of the strong 
dependence of A^o on that splitting Q . 

Results from Deltuva and Fonseca using the Doleschall 
potential INOY04 j4l| are also shown in Figures!^ andlTUl 
That potential introduces non-localities to simulate im- 
plicitly the effect of three-nucleon forces, which are nec- 
essary to reproduce three-and four-nucleon binding ener- 
gies. The parameters of those non-localities are adjusted 
to reproduce 3N scattering phase-shifts, and also better 
to reproduce the AyQ measurements of Ref. . Here, this 
potential improves the description of Aqj, somewhat, but 
has little effect on A^x and Ayy. Considering the phase- 
shifts, this model better describes ^Pq and produces a 
P-wave splitting closer to the experimental results. 

The addition of explicit phenomenological 3N forces 
in theoretical calculations has traditionally not provided 
full agreement between experiment and theory for p -|- 
•^He observables, especially for AyQ. New results from 
Viviani [H, HI] using a 2N [131 and 3N Q interaction 
derived from chiral perturbation theory at N3L0 and 
N2L0, respectively, are shown in Fig. [TU] The calcula- 
tions, made using the Kohn variational principle and the 
hyperspherical harmonic technique, are compared with 
both the present Aq^ and Aj^j^ results and Aj,o from Ref. 
[3| at 4 MeV. Satisfying reduction of Aj^ puzzle differences 
was obtained using this effective-field-theory version of 
the 3N interaction, corresponding to the improved agree- 
ment with the triplet P-wave phase-shifts shown in Fig. 
ini Better agreement for the e(l~) mixing parameter is 
also evident. Though theoretical agreement with exper- 
imental results is still not complete, Machleidt has sug- 
gested that sizeable one-loop 3N force diagrams exist 
at N4L0 of the A-less chiral theory, or at N3L0 when a 
phenomenological A is included, and that their addition 
may ultimately explain the remaining differences. 



V. CONCLUSIONS 

In this paper we have presented new measurements 
of Aoy, Axx, and Aj^j^ for p-l-'^He elastic scattering be- 
tween 2 and 6 MeV proton energy. The target analyz- 
ing power measurements represent an improvement in 
accuracy over previous results, while the spin-correlation 
measurements include the lowest-energy data to date. 

These new measurements were included in a global 
phase-shift analysis and new phase-shifts were extracted. 
These additional data remove the ambiguity reported 
in Ref. @. Though discontinuities in the energy- 
dependence of the D- and F-wave phase-shifts were 
present, a single global solution was obtained by requir- 
ing that all phases be continuous over the energy range 
of the global database. 

Recent theoretical calculations [s^l using realistic 2N 
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FIG. 9: (Color online) Comparison of the present p-wave 
phase-shifts, in degrees, with the theoretical results of Deltuva 
and Fonseca [37] using the CD-Bonn and Doleschall INOY04 
potentials and of Viviani ;43| using a potential derived from 
chiral potential theory (ChPT). The splitting S is defined in 
the text. 



potentials underpredict the present Aq^ results by 10- 
20%, which is similar to but smaller than the previously 
observed 40% underprediction of A^o [1, 113] • The spin- 
correlation coefhcients A^x and Ayy are better-described, 
though small underpredictions are observed at forward 
angles, and qualitatively different trends are observed at 
2.25 MeV. The S-wave phase shifts agree well, but ^Pi, 
■^Pq, and e (1^) differ, and the theoretical triplet P-wave 
splitting is too small. 

The INOY04 potential [4l| which includes a phe- 
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FIG. 10: (Color online) Comparison of the present observ- 
ables Aoj/ and Ayy measured at 4 MeV, along with the A^o 
measurements of Fisher et al. [j, with the theoretical re- 
sults of Viviani [i^] using a using a potential derived from 
chiral potential theory (ChPT). The results of Deltuva and 
Fonseca [sj using the CD-Bonn and INOY04 potentials are 
also shown. 



nomenological 3N force improves the description of Aqj, 
and '^Po and increases the P-wave splitting. Preliminary 
calculations by Viviani [42] of p + '^He elastic scatter- 
ing observables at low energies using chiral 2N and 3N 
potentials show closer agreement with A^o- 
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